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We show how the appearance of d-wave pairing in fermionic condensates manifests itself in inelastic 
hght scattering. SpecificaUy, we calculate the Bragg scattering intensity from the dynamic structure 
factor and the spin susceptibility, which can be inferred from spin flip Raman transitions. This 
information provides a precise tool with which we can identify nontrivial correlations in the state 
of the system beyond the information contained in the density profile imaging alone. Due to the 
lack of Coulomb effects in neutral superfluids, this is also an opportunity to observe the Anderson- 
Bogoliubov collective mode. 

PACS numbers: 03.75.Fi,67.85.De,74.25.-q,03.75.Nt 
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INTRODUCTION 

Ultra-cold atomic systems provide a way to build quan- 
tum simulators - simulating a model with an experi- 
mental system [l|. Studying models for strongly corre- 
lated systems can shed light on outstanding problems 
in condensed matter physics, for example whether the 
ground state of the doped Hubbard or t-J models sup- 
ports d-wave superconductivity [2]. It is important to 
have precise tools to identify and characterize the result- 
ing phases. Establishing unconventional pairing symme- 
tries requires new techniques in the setting of ultra-cold 
gases. Traditional probes like phase sensitive measure- 
ments, transport, or scanning tunneling microscopy, are 
all unavailable to condensates in optical systems. How- 
ever, techniques do exist which provide information be- 
yond the scope of what time-of- flight imaging can do, 
allowing precise identification of new phases. 

We calculate the zero temperature response of a 
fermionic d-wave superfiuid to inelastic light scattering. 
The primary emphasis in our work is that the inelas- 
tic response functions we provide are specific enough to 
unambiguously identify the underlying d-wave symme- 
try. We also have the opportunity to observe collective 
modes in a neutral paired superfiuid. To this end, we 
also look at s-wave pairing. In solids, this mode is lifted 
to the plasma frequency by the Coulomb interaction and 
disorder can play an important role, making the collec- 
tive mode difficult to observe. In the context of conden- 
sates, the complications from the Coulomb interaction 
and disorder are absent. A disadvantage is the effect of 
the confining potential which is usually present even for 
optical lattices. 

Previous studies have focused on identifying s-wave su- 
perfiuidity [3|47|, or using noise-correlations to identify 
many-body states [8] , or using a periodic driving of the 
system then observing its response [9] . Collective modes 
in chiral p-wave superconductors have also been exam- 
ined [10]. We explore the confiuence of these ideas, fo- 



cusing on identifying d-wave superfiuids, and on the effect 
of their collective mode. 



BRAGG SCATTERING 



We recall the textbook [ll| discussion of inelastic light 
probes here for clarity. Bragg scattering is a two photon 
process that transfers momentum Q = ki— k2 and energy 
Q = uji —UJ2 to the ground state, and can be used to probe 
density ffuctuations of the system. The probing light can 
be tailored to couple to density and between hyperfine 
degrees of freedom. The intensity of the light from the 
laser beams takes the form, /q cos(Q • r — Qt)^ and the 
atoms in the sample will feel an optical potential due 
to the AC Stark effect. We label the couplings to the 
density (d) or spin (s) degree of freedom as Vd^s- The 



perturbation is taken as JjJ, ll2| 
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where Spd^s,Q is the variation in the density or spin den- 
sity. We will generalize this later in the paper. 

The dynamic structure function 6'(Q, Q) is the Fourier 
transform of the density-density correlation function, 
where Q and Q are respectively the momentum and 
energy transferred by the probe. The excitation rate 
per particle is given by 27r{^)'^S{Q,n) [l2J. The dy- 
namic structure factor can be measured in more than one 
way 12|, ll3|. For example, after the atomic sample is il- 
luminated by the two laser beams, the trap is switched 
off, and time of ffight images can measure the Bragg scat- 
tered atoms - distinguished because they are displaced 
by the absorbed momentum. What is actually mea- 
sured in this experiment is the rate of momentum transfer 



which is ^ = Q{^)^^[S{Q,n) - S{-Q,-n)] ex 
3-^^)^(Q, n) -^ \imT=o S{Q, n). In the end, Bragg 
scattering provides a measurement of the imaginary part 
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of the density-density correlation function which we cal- 
culate. Similarly, it has been suggested that the spin flip 
rate observed by Raman transitions can be used to study 
the spin susceptibility [6|, ll3| . These two probes are spe- 
ciflc enough to conflrm unconventional Fermi superfluid 
pairing when it arises in condensates. 



HAMILTONIAN AND RESPONSE FUNCTIONS 

We limit our discussion to a dilute Fermi gas which un- 
dergoes pairing in the BCS limit. The BCS Hamiltonian 
generalized for anisotropic pairing is: 



^ = Y1 (^kt ^-j^^) 
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(2) 



For a two dimensional square lattice in the tight bind- 
ing approximation, the dispersion is 



^k = -2t[cos{kx) + cos{ky)] 
The d-wave superconducting gap is 



/i. 



(3) 



Ak = AoXk, Xk = cos{ka:) - cos{ky). (4) 

We chose the chemical potential to be ja = — t, and 
Aq = 0.2t. The resulting Fermi surface is shown in Fig.[TJ 
We note that particle-hole symmetry takes fi = —t ^ 
ji = t and is just a translation of the Fermi surface by 
(7r,7r). The Bogoliubov quasiparticles have an excitation 
spectrum given by: E\^ = y^S,^. + A^. The Green's func- 
tion for a d-wave superconductor takes the well known 
form in Nambu space. 



G(k, iUn) = (icJni + ^kTs + Akfi) = 
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FIG. 1: The Fermi surface for a square lattice in the tight 
binding limit including only nearest neighbor hoppings at /i = 
—t. The nodal lines and sign of the gap are marked. 
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FIG. 2: (Color Online) Plot of Im[n22(0, 0) - n22(Q, ^)]"^ 
in the s-wave state for (a) Q = 0.5Ao (b) Q = 1.5Ao (c) 
Q = 2.OA0. The poles in the denominator of the second term 
in Eq. (|18p illustrate the Anderson-Bogoliubov mode. Panel 
(b) shows the approximate dispersion captured from panel 
(a) in the (1,0) and (1,1) directions compared to the result 



In this work, the inclusion of both the d-wave quasipar- 
ticle and collective response is crucial. Following the pio- 
neering works of Kulik et al. Il5|, Wong and Takada [l6|, 
and Ohashi and Takada [l7|, ll8| , we calculate the phys- 
ical response by including fictitious interactions. To ac- 
complish this we generalize the density operators. The 
Pauli matrices in Nambu space will be denoted by fi with 
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c_k|) is a spinor in Nambu space. 



In this notation, the generalized densities are 



PiQ 



X]7k+Q/2, 
k 



^If,^k+Q, 



7k,i=0, 1,2,3 = (l,Xk,Xk,l), 



P^Q = Y^ 4,t^k+Q,t - Ck,i^k+Q,i (6) 

k 

PlQ = Y. ^k+Q/2 (4t^-k-Qi + h'^) (7) 

k 

P2Q = Y ^k+Q/2 (^4t^-k-Qi + ^-C-) (8) 

k 

ha = Yl 4,t^k+Q,t + 4,i^k+Q,i. (9) 

k 

Physically, psQ = pdQ is the density operator, while 
piQ and P2Q can be thought of as, respectively, qual- 
itatively encapsulating the amplitude and phase of the 
order parameter. The spin-density, treated separately 
throughout the paper, is psQ = poQ- By treating the 
spin density separately we do not include the possibility 
of spin-phase, spin-amplitude, etc. modes. According to 
Buchler et al. [6], the contribution from collective terms 
for /TTipoo] vanish, so /TTipoo] = /mpQQ]. In Zou et 
al.[19| a self-consistent spin-spin interaction is included 
in the Random Phase Approximation for the s-wave case 
which goes beyond our analysis. For interactions, we in- 
clude the terms: 



Hint = ^(PlQPl,- 



-Q +P2QP2,-q). 



(10) 



The density and spin response of the system is calcu- 
lated using linear response theory. The density-density 
response has been generalized to be a matrix of den- 
sity responses in Nambu space. The generalized response 
functions, 






take the form: 
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The gap equation at zero temperature, 
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is equivalent to = 1 + f n^2(0. 0) where "g" is the BCS 
coupling constant stemming from assuming a separable 
pairing interaction F(k, k') = ^XkXk'- Zero temperature 
Bragg scattering would correspond to a measurement of 
Imn33(Q,l]), and a different experiment could measure 
Imnoo(Q,^), the spin susceptibility. 

The density-density response at T=0 comes from the 
following integrals 

— ^da;7k+Q/2,i7k+Q/2,i 



uj{uj + 1^) + a^k^k+Q + 6AkAk+Q 
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(14) 



The upper indices are a, 6 = ±1 and are determined by 
which Nambu-space matrices occur in the lower indices 
for the generalized densities. The spin susceptibility {i = 
0^=0) corresponds to a = 1, b = 1 and has been 
previously calculated [20]. The density channel (i = 3 
j = 3) has a = 1, b = — 1 and the phase channel {i = 2 
j = 2) has a= — l,b = —1. After completing the contour 
integral, it is helpful to replace the terms E{E -\- Q) by 
a choice which cancels some terms in the denominators, 
E{E ^n) = E{E ^E' ^n)- EE'. 
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(15) 



Q + £^k+Q -^ E\^^ir] £^k+Q + £^k - ^ + ir]J 

The phase-density response function 1132 = —1123 is 
also needed: 

rr ^o o^ o- f _f^ Ak£^k+Q + £^kAk+Q^ ^^^ 
n32(Q,^) = 2.y (^7k+Q, ^^^^^^^ (16) 

1 1 



ft + £^k+Q -^ Ei^^if] £^k+Q -^ Ei^-ft^if] 



The Feynman diagrams of the random phase approxi- 
mation (RPA) are resummed for a neutral superfluid. It 
is found by explicit calculation |15l4l7| that the amplitude 
modes decouple from the density and phase oscillations, 
ni2 = ni3 = 0, so the the density response is: 



( 



n22 n23 
\ 1132 1133 



^-^22 ^-^23 
^^32 ^^33 



-g/2 




The full density-density response [16], after analytic 
continuation Qm ^ ^ -\- if] (for numerics r] was 10~^t), 

Note that the gap equation is used to eliminate g in 

termsof n^2(0.0)- 

The formula in Eq. (p!8|) has a clear interpretation. The 
first term is the ordinary density-density response. The 
second term is a coupling between the density and phase 
response with a phase mode in the denominator. A 
pole in the second term is the Goldstone, or Anderson- 
Bogoliubov, mode for the neutral superfluid. 

The finite temperature expression is required to obtain 
the normal state limit. This T=0 result only contains the 
pair-breaking terms, because the scattering terms vanish 
for T=0, and both are necessary to recover the normal 
state limit, which follows directly from the result |2l| . 

The calculation of the inelastic neutron response, no- 
tably in cuprates [22|, |23| , is very similar to this calcu- 
lation; frequently, in cuprate calculations J2J, |25| the 
RPA analysis contains a magnetic structure J(q) = 
[cos{qx) + cos{qy)]^ bilayer effects, or Coulomb interac- 
tions which are absent here. 



DISCUSSION 

If the system were to condense into an s-wave super- 
fluid, we should see the effect of opening a gap on the 
Fermi surface. The only low-energy feature present in 
the s-wave case is the Anderson-Bogoliubov mode. All 
non-Goldstone excitations are gapped below the scale 
Q < 2Ao. In Fig. [2^ we show the remaining gapless 
collective excitations. Maxima in the signal for various 
frequencies are compared to the exact result Q = ^Q 
[26| in Fig.[2]3. The presence of a square lattice dispersion 
causes anisotropy in the Anderson-Bogoliubov mode's 
dispersion near the pair breaking frequency Q = 2Ao 
shown in Fig. [2j In ordinary s-wave superconductors, 
the presence of the Coulomb interaction is said to lift the 
energy of this excitation to the plasma scale. In a neu- 
tral superfluid, the absence of a Coulomb interaction and 
gapped quasiparticle excitations means that the only low 
energy mode is this collective mode. 

The response functions calculated in this paper are 
compHcated "speckle" patterns in Q-space, shown in 
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FIG. 3: (Color Online) The imaginary part of the zero temperature spin susceptibility, Impoo], Eq.[T2]in the d-wave state, for 
frequencies [a)Q. — 0.5Ao, (&)Ao, (c)1.5Ao, with Ao = 0.2t. Notice that there is zero response at Q = (0,0) consistent with 
the fact that the spin structure factor vanishes in the zero momentum limit [y]. 
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FIG. 4: (Color Online) The imaginary part of the zero temperature density-density response, Impss], Eq. [T8]in the d-wave 
states, for frequencies (a)Q = 0.5Ao, (^)Ao, (c)1.5Ao. We have artificially put a mark in the center of Figs (a) and (c) where 
the data was clipped to represent points off the scale. For these plots Ao = 0.2t. 



Figs. [3] and m Define the coherence factor 77^5 (k,Q) 
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, Kj ^ ^ ,. The response function can 

^ ^k^k+Q ^k^k+oy ^ 

be interpreted from Eq. [15] as a convolution proportional 
to: 

Im n^^ (X / duiduj2 Y^ ^ab(k, Q) 



X d{Q -uji- uj2)S{uJi - E\^)S{uj2 - ^k+q). (19) 

The response function is a coherence-factor-weighted 
sum over the delta function. The single particle den- 
sity of states (DOS) is DOS(a;) = -^Y.k^^G{k,uj) = 
"^^k ^(^ ~ -^k)- If the frequency integrals in Eq. [19] are 
performed then we can also define a joint density of states 
(JDOS), for two particle excitations as 

JDOS(Q, n) = Yl ^(^ - Ek+Q - Ek). (20) 



By examining the structure of the momentum resolved 
DOS and JDOS, meaning DOS(k,cj) = S{uj - Ek) for 



fixed CO without a sum over k, and using reasoning anal- 
ogous to the octet model [271, |28| from cuprate tunnel- 
ing experiments [29| , the structure of the response func- 
tions becomes clear. We show plots of DOS(k, a;) for a 
d-wave superconductor defined by Eq.(j5|) in Fig. [5] The 
essential observation is that the final observed signal at 
wave vector Q will have a high transition probability 
when DOS(k,cj) = DOS(k + Q,cj + Q). So we com- 
pare the two densities of states. The available states in 
the DOS(k, uj) convolve to give rise to the JDOS (Fig. [6]), 
which is then, in turn, weighted by a coherence factor in 
the dominant contribution to the response function. In 
tunneling experiments the scattering was elastic in con- 
trast to this inelastic process transferring energy Q from 
the state at E\^ to the state at £^k+Q, so we compare 
DOS(k,cj) = DOS(k + Q,uj). This differs from the in- 
elastic case which blurs the DOS over a range of allowed 
frequencies. 

In Fig. [6] we have overlayed the JDOS for several 
fixed external frequencies with the response function. 
There is an overall qualitative agreement between the Q 
which nest the single DOS(k, u) and DOS(k + Q, cj + 1]) 
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FIG. 5: (Color Online) Momentum resolved density of states 
for fixed frequency Q. For Q = O.lAo(a) we see nodal excita- 
tions. For Q = Ao (b) the region of highest DOS corresponds 
minima in Vk^k- The solid black line is the equal energy 
contour for the same energies. 



and the JDOS's structure. The gross features of the 
JDOS (Fig. [6]) are contained in the final response func- 
tion (Figs. [3] and H]) which differ by a coherence factor. 
In both figures [3] and [3J we see nodal excitations low 
(1], Q) along the diagonals as well as transitions at finite 
Q resulting from the convolutions of the DOS. 

The collective mode, shown in Fig. [3 has a structure 
beyond the simple s-wave result due to it being d-wave. 
To check this, the poles in second term in Eq. ([18]), 
Im[n22(0,0) - n22(Q,^)]"\ can be plotted by itself 
in Fig. [71 These excitations are the d-wave equivalent 
of the Anderson-Bogoliubov mode. Unlike the s-wave 
Anderson-Bogoliubov mode where all quasiparticle exci- 
tations are gap, the d-wave Anderson-Bogoliubov mode 
is not isolated from the gapless quasiparticle response in 
the full Impaa]. 

An additional observation can be made. In the in- 
elastic neutron response for cuprates, excitations around 
(7r,7r) have played an important role ^]. We chose not 
to focus on the (0,7r) excitations which are broader and 
weaker than the (tt, tt) excitations. Consequently, we are 
curious to see what comes out of the RPA analysis per- 
formed here. We have tracked the peak excitations which 
occur around the point (tt, tt) for a limited range of fre- 
quencies. Within the limitations imposed by sampling 
discrete frequencies, we find that clear peak excitations 
around (tt, tt) are absent for low energy, begin along the 
zone diagonal, and split into two peaks as we increase 
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FIG. 6: (Color Online) Joint density of states, JDOS(Q,^), 
for two-particle excitations over the range (0,7r) (a,b) and con- 
tour plots of the same (c,d). 



frequency. At the pair-breaking energy 2Ao = 0.4t those 
excitations are near the Brillouin zone boundary, and 
above that energy they are difficult to identify amidst all 
the allowed excitations. This is shown in Fig. [H 

We have shown by taking apart Eq.[T9]that a combina- 
tion of the geometry of the Fermi surface and the quasi- 
particle dispersion determines the "speckle" patterns in 
Figs. [3] and [H Notwithstanding the complicated na- 
ture of both the spin susceptibility and density-density 
response, the detailed k-space features should be specific 
enough to identify exotic pairing in superfluids distinctly 
from other states. 



CONCLUSION 

In summary, we have provided an analysis of two re- 
sponse functions, the spin susceptibility and the density- 
density response, for the inelastic scattering of light off a 
d-wave condensate. 

We included the effect of collective modes and exam- 
ined their dispersion in the s-wave and d-wave cases. In 
the s-wave case, the collective mode is the only mode 
present below the pair breaking energy 2Ao, unencum- 
bered by effects from the Coulomb interaction and dis- 
order in the context of cold atom experiments. In the 
d-wave case, the collective mode is superimposed with 
other gapless excitations and has a structure beyond the 
simple s-wave case. 

For the d-wave spin and density reponse, we provide 
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FIG. 7: (Color Online) Plot of /m[n22(0,0) - n22(Q,r2)]"^ in the d-wave state for (a) Q = Ao (b) Q = 1.5Ao (c) Q = 2.0Ao. 
The poles in the denominator of the second term in Eq. (|18p illustrate the Anderson-Bogoliubov mode. This result differs from 
the s-wave Anderson-Bogoliubov mode, Fig. [2l due to the presence of a d-wave gap. 
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FIG. 8: (Color Online) Positions of the maxima in Imllss 
around the point (tt, tt) in Q — Q space. 



an interpretation of the detailed momentum-space struc- 
ture. The inelastic response functions reflect the opening 
of a d-wave gap on top of the geometry of the Fermi 
surface. These patterns are specific enough to unam- 
biguously identify d-wave superfluidity, were it to arise, 
providing a cold-atom analog to phase-sensitive exper- 
iments in solids. We also examined the peak response 
near the point Q = (tt, tt) in analogy to the hourglass in 
cuprate neutron experiments [30] . The primary result of 
our work is a precise tool, the momentum space structure 
of inelastic light scattering, to identify and characterize a 
d-wave paired Fermi superfluid which can be generalized 
to any gap symmetry. 
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